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ABSTRACT 

Discoveries from the prime Kepler mission demonstrated that small planets (< 3 I?®) are common 
outcomes of planet formation. While Kepler detected many such planets, all but a handful orbit 
faint, distant stars and are not amenable to precise follow up measurements. Here, we report the 
discovery of two small planets transiting K2-21, a bright (K = 9.4) MO dwarf located 65 ± 6 pc from 
Earth. We detected the transiting planets in photometry collected during Campaign 3 of NASA’s 
K2 mission. Analysis of transit light curves reveals that the planets have small radii compared to 
their host star, Rp/R * = 2.60 ± 0.14% and 3.15 ± 0.20%, respectively. We obtained follow up NIR 
spectroscopy of K2-21 to constrain host star properties, which imply planet sizes of 1.59 ± 0.43 J?® 
and 1.92 ± 0.53 I?®, respectively, straddling the boundary between high-density, rocky planets and 
low-density planets with thick gaseous envelopes. The planets have orbital periods of 9.32414 days 
and 15.50120 days, respectively, and a period ratio P c /Pb = 1-6624, very near to the 5:3 mean motion 
resonance, which may be a record of the system’s formation history. Transit timing variations (TTVs) 
due to gravitational interactions between the planets may be detectable using ground-based telescopes. 
Finally, this system offers a convenient laboratory for studying the bulk composition and atmospheric 
properties of small planets with low equilibrium temperatures. 

Subject headings: K2-21— EPIC-206011691— techniques: photometric —techniques: spectroscopic 


1. INTRODUCTION 

Analysis of photometry collected during the prime Ke¬ 
pler mission (2009-2013) demonstrated th at small plan¬ 
ets are common around G , K, fc M dwarfs dHoward et al. 
2012 IFressin et al]|2013t IPetigura et al.||2013a| [Dress¬ 
ing & Charbonneau 2013a L Planets that are nearly 
Earth-size (1-2 A®) and orbit close to their host stars 
(P < 100 days) are abundant around Sun-li ke stars (26% 
of G fe K dwarfs host such a planet; Petigura et al.| 
2013a) and near ly ubiquitious around M dwar fs (1.2 
planets per star; Dressing & Charbonneau 2015). The 
higher occurrence of such planets around M dwarf hosts 
cannot be explained by higher detection efficiency for 
small planets around M dwarfs. Both Petigura et al.| 
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(2013a) and Dressing & Charbonneau (2015) empirically 
derived survey completeness by injecting mock transits 
into Kepler photometry and measuring the recovery rate. 
Both groups corrected for modest pipeline incomplete¬ 
ness (< 50% for P < 100 days, Rp = 1-2 R®). An open 
question is whether small planets are intrinsically more 
numerous around M dwarfs, or if M dwarf systems are 
more compact. 

M dwarfs offer a convenient laboratory to study the 
bulk physical properties and atmospheres of small plan¬ 
ets. Planet transits are deeper and radial velocity sig¬ 
natures are larger, for a given planet size and orbital 
period. Planets around M dwarfs having low equilib¬ 
rium temperatures (e.g. T eq = 200-600 K) have tighter 
orbits with shorter orbital periods compared to simi¬ 
lar planets around solar-type stars. These planets are 
especially compelling targets for atmospheric transmis¬ 
sion spectroscopy by the James Webb Space Telescope 
(JWST), provided a sample of nearby, bright M dwarfs 
with warm transiti ng planets in the ~ l-3 i?® size range 
can be identified (Batalha et al. 2015). 

Finding planets that transit M dwarfs is difficult be¬ 
cause the observable stars are generally faint, and the 
brightest (nearest) ones are sparsely distributed on the 
sky. Kepler surveyed a few thousand M dwarfs and 
only discovered «160 planet candidates, despite the 
high planet occurrence rate. To capitalize on the fol¬ 
low up opportunities, bright stars are required. Al¬ 
though the Kepler sample of M dwarfs was faint (me¬ 
dian Kp = 15.5 mag), the magnitude-limited bright¬ 
ness distribution naturally included some brighter tar¬ 
gets. Among the most favorable is Kepler-138 (Kp = 12.9 
mag, K = 9.5 mag), a system of three transiting planets, 
includi ng two nearly Earth-size planets (1.2 R m ) in wa rm 
orbits ( Jontof-Hutter et al.||2015 Kipping et al. [2014 1. 
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Now that Kepler is operating in the K2 mode, we have 
the opportunity to greatly expand the number of known 
small planets transiting bright M dwarfs. During K2 ob¬ 
servations, the Kepler Space Telescope observes different 
regions of the ecliptic every 90 day s, castin g a wider net 
for nearby transiting planets (Howell et al. 2014) . Essen¬ 
tially, each of the planned fourteen K2 pointings offers 
the possibility of discovering a sample of short-period 
planets orbiting bright stars that is similar to the ensem¬ 
ble from the prime Kepler mission. 

Here we present the discovery of a new multi-planet 
system orbiting a bright M dwarf in the K2 Campaign 3 
field, K2-21, also known as EPIC-206011691. This sys¬ 
tem is a cousin of Kepler-138, identified during the prime 
Kepler mission, and K2-3, a system of three transiting 
super-Eart hs discovered by our t eam in an earlier K2 
Campaign (Crossfield et al. 2015). It also offers a pre¬ 
view of expected results from the TESS mission, when 
sky cover a ge will i ncrease by another order of magnitude 
(Ricker et al. ||2015 1. 

In this paper, we describe our detection of K2-21b and 
c from K2 photometry along with our spectroscopic and 
imaging follow up in Section [2] In Section [3j we present 
our analysis of stellar properties, planet properties, and 
false positive probability. Finally, we place the K2-21 
system in the context of other transiting planets in Sec¬ 
tion [4] and offer some thoughts regarding the rich follow 
up potential for this system. 

2. OBSERVATIONS 

2.1. K2 Photometry 

K2-21 was observed during K2 Campaign 3 lasting 
from 2014 November 14 to 2015 February 3. The star 
is listed as EPIC-206011691 in the Mikulski Archive for 
Space Telescopes (MAST). Target proper ties, inclu ding 
optic al and NIR pho tometry from APASS (Henden et al. 
20121) and 2MASS (Skrutskie et al. 20061 are listed in 
Table [l] 

We extracted K2 photometry for the star K2-21 from 
the Kepler pixel data, which we downloaded from the 
MA ST. Our photometric extraction routine is outlined in 
Crossfield et al. (2015). In brief, during K2 observations, 
stars drift across the CCD by ~1 pixel every ^6 hours. 
As the stars drift through pixel-phase, intra-pixel sen¬ 
sitivity variations and errors in the flatfield cause the ap¬ 
parent brightness of the target star to change. For every 
observation we solve for the roll angle between the target 
frame and an arbitrary reference observation. We model 
the time- and roll-dependent brightness variations using 
a Gaussian process. Figure l] shows both the raw and 
corrected photometry. We do not see any evidence for 
stellar variability with timescales of < 5 days, suggest¬ 
ing that K2-21 is a slowly rotating star. The corrected 
light curve exhibits ~1% variability on timescales longer 
than ~10 days. This provides an upper limit on the 
intrinsic variability of K2-21 over long timescales. How¬ 
ever, since stars drift perpendicular to the roll direction 
over the course of a campaign, it is difficult to disen¬ 
tangle long-term astrophysical variability from position- 
dependent variability. Our calibrated K2 photometry is 
included as an electronic supplement. 

2.2. Transit Detection 


We searched through the calibrated and detrended 
photometry (top pane l of Figure [2]) using th e TERRA 
algorithm described in Petigura et al. (2013b I. TERRA 
identified a transit candidate having P = 15.50120 days 
and signal-to- noise ratio (SNR) = 25. We fit this can¬ 
didate with a Mandel & Agol (2002) model and sub¬ 
tracted the best fit model from the photometry. We reran 
TERRA on the photometry with the P = 15.50120 day 
candidate removed. We found a second candidate having 
P = 9.32414 days and SNR = 22. Again we removed the 
best-fitting model. In subsequent searches, TERRA did 
not find any additional periodic box-shaped signals. 

2.3. Infrared and Optical Spectroscopy 

We observed K2-21 on 2015 July 23 UT using the 
recen tly refurbished SpeX spectrograph (Rayner et al. 
2003) on the 3.0 m NASA Infrared Telescope Facility 
(1K1F). The data were taken under clear skies with an 
average K-band seeing of 0.4-0.5 arcsec. We observed 
with the instrument in short cross dispersed mode (SXD) 
using the 0.3 x 15 arcsec slit. This setup provides simul¬ 
taneous wavelength coverage from 0.7 to 2.5 /mi at a 
resolution of R ss 2000. The extended blue wavelength 
coverage is a result of the recent chip upgrade SpeX re¬ 
ceived in 2014 July. The target was placed at two po¬ 
sitions along the slit and observed in an ABBA pattern 
for subsequent sky subtraction. The observing sequence 
consisted of 8 x 60 s exposures for a total integration time 
of 480 s. Once the exposures were stacked, this integra¬ 
tion time led to a signal-to-noise of ~ 200 per resolution 
element. We obtained standard SpeX calibration frames 
consisting of flats and arclamp exposures immediately 
after observing K2-21. 

The S peX spectrum was reduced using the Sp eXTool 
package (IVacca et al.||2003| |Cushing et al.||2004). SpeX¬ 
Tool performs flat fielding) bad pixel removal, wave¬ 
length calibration, sky subtraction, spectral extraction 
and combination, telluric correction, flux calibration, and 
order merging. Flux calibration and telluric corrections 
were perfomed using the spectrum of the AOV-type star 
HD 218639 which was observed within 16 minutes and 
0.04 airmass of the target. The calibrated spectrum is 
compared to late-type standards from the IRTF Sp ec- 
tral Library (Cushing et al.||2005l Rayner et al.||2009 1 in 
Figure [3j K2-21 is a best visual match to the K7/M0 
standaras across the near-IR band s. W e perform a more 
detailed spectral typing in Section |3.1| 


K2- 21 was also observed using Keck/HIRES (Vogt 


et al. 1994) on 2015 July 24 UT with a total integra¬ 
tion time of 216s. HIRES provides wavelength cov¬ 
erage from ~3600-8000 A at R ze 60000. The spec- 


trum was reduced using the standard 

pipeline of the 

California Planet Search (Marcy et al. 

2008 

1 and has 


2.4. Archival and Adaptive Optics Imaging 

To assess the possibility of background stars falling in¬ 
side the software aperture of Kepler pixels, we compare 
two epochs of imaging data from the Digitized Sky Sur¬ 
vey (DSS). In Figure E] we show the DSS-Blue plates 
taken on 1954 July 25 (top panel) and the DSS-Red 
plates taken on 1991 July 17 (bottom panel). The images 
are 2x2 arcmin and have a pixel scale of 1 arcsec/pixel. 
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TABLE 1 

Stellar Parameters of K2-21 


Parameter Units Value Source 


Identifying information 


EPIC ID 

- 

206011691 

EPIC 


2MASS ID 

- 

22411288-1429202 

2MASS 


a R.A. 

h:m:s 

22:41:12.89 

EPIC 


S Dec. 

d:m:s 

-14:29:20.35 

EPIC 


i 

d:m:s 

48:56:11.39 

EPIC 


b 

d:m:s 

-57:10:44.78 

EPIC 


Photometric Properties 



Kp 

mag 

12.31 

EPIC 


B 

mag 

14.14 ± 0.06 

APASS 


V 

mag 

12.85 ± 0.02 

APASS 


g' 

mag 

13.53 ± 0.02 

APASS 


r' 

mag 

12.32 ± 0.06 

APASS 


i' 

mag 

11.76 ± 0.06 

APASS 


J 

mag 

10.25 ± 0.02 

2MASS 


H 

mag 

9.63 ± 0.02 

2MASS 


Ks 

mag 

9.42 ± 0.02 

2MASS 


Spectroscopic and Derived Properties 



f-f'Ot 

mas yr —1 

17.3 ± 1.4 

Zacharias et al. 

(2012 

AM 

mas yr —1 

-78 ± 2.8 

Zacharias et al. 

(2012 

Spectral Type 

- 

M0.0T0.5 

SpeX, this paper 

T e ff 

K 

4043 ± 375 

SpeX, this paper 

[Fe/H] 

dex 

-0.11 ± 0.13 

SpeX, this paper 

M* 

M© 

0.64 ±0.11 

SpeX, this paper 

it* 

Rq 

0.60 ± 0.10 

SpeX, this paper 

P* 

g cm -3 

4.15 ±0.61 

SpeX, this paper 

L* 

L© 

0.086 ± 0.064 

SpeX, this paper 

Distance 

pc 

65 ±6 

this paper 


Age 

Gyr 

>1 

this paper 




Fig. 1.— Top : Raw photometry computed by summing the background-subtracted counts inside a circular aperture (3 pixel radius) 
centered on K2-21. Bottom : Photometry after correcting for variations due to telescope roll angle. Noise on three-hour timescales has been 
reduced by a factor of 30. The ~1% variability gives an upper limit to K2-21’s intrinsic variability. However, since stars drift perpendicular 
to the roll direction over the course of a campaign, it is difficult to disentangle long-term astrophysical variability from position-dependent 
variability. The data used to produce the bottom panel is included as an electronic supplement. 

The images are centered on the epoch 2015 coordinates of the target as observed by Kepler (a: 22:41:12.9, <5: 
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Fig. 2.— Calibrated K2 photometry for K2-21. Vertical ticks indicate the locations of each planet’s transits. Bottom: Phase-folded 
photometry and best-fit light curves for each planet. Best fit parameters from light curve fitting are tabulated in Table [2] 
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Fig. 3;_ Calibrated IRTF/SpeX spectra of K2-21 compared to spectral standards. After fitting several spectral indices (described in 

Section ED- we derive a best-fit spectral type of M0.0d=0.5. Spectroscopically derived stellar parameters are listed in Table [I] 


-14:29:21.9 J2000.0) and the open circle represents the 
aperture size used when extracting the calibrated pho¬ 
tometry. 

The DSS images clearly show the proper motion of 
the primary target, while the nearby stars, located 40" 
W and 60" SE, show no significant astrometric motion. 
The primary target, in contrast, displays a clear proper 
motion of ^3" over 37 years, in reasonable agreement 


with the measured prop er motion (Lepine & Gaidos 2011 


Zacharias et al. 2012). In the JJ55 image there is no 
evidence ot a background star at the 2015 position of K2- 
21. We estimate that if a star is located this position, 
that star must be at least 2.5-3 mag fainter than K2-21. 

We also obtained near-infrared adaptive optics images 
of EPIC 206011691 at Keck Observatory on the night of 
2015 July 24 UT. Observations were obtained with the 
1024 x 1024 NIRC2 array and the natural guide star 
system; the target star was bright enough to be used as 
the guide star. The data were acquired using the Kcont 
filter using the narrow camera field of view with a pixel 
scale of 9.942 mas/pixel. The Kcont filter has a narrower 


bandwidth (2.25-2.32 pm) compared the K filter (1.99- 
2.40 /im) and allows for longer integration times before 
saturation. A 3-point dither pattern was utilized to avoid 
the noisier lower left quadrant of the NIRC2 array. The 
3-point dither pattern was observed thee times with 10 
coadds and a 1.5 second integration time for a total on- 
source exposure time of 3 x 3 x 10 x 1.5s = 135 s. 

The target star was measured with a resolution of 
0.06 arcsec (FWHM). No other stars were detected 
within the 10 arcsec field of view of the camera. In the 
Kcont filter, the data are sensitive to stars that have K- 
band contrast of A A' = 3.7 at a separation of 0.1 arcsec 
and AK = 7.5 at 0.5 arcsec from the central star. We 
estimate the sensitivities by injecting fake sources with 
a signal-to-noise ratio of 5 into the final combined im¬ 
ages at distances of JVx FWHM from the central source, 
where N is an integer. The 5<r sensitivities, as a function 
of radius from the star, are shown in Figure [5j 

3. ANALYSIS 
3.1. Stellar Parameters 
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Fig. 4.— POSS1 blue plates observed in 1954 (top panel) and 
POSS2 red plates observed in 1991 (bottom panel). The circle 
shows the location and extent of the circular aperture used to ex¬ 
tract K2-21 photometry at the 2015 position of the star. Between 
1954 and 1991, K2-21 moved by ~3 arcsec, which can be clearly 
seen in the DSS images. In the POSS1 plate, the star is offset from 
the 2015 position by ~ 5 arcsec. The DSS blue plates rule out a 
background star coincident with the current location of K2-21 to 
A B = 2.5-3.0. 

Here, we derive K2-21 effective temperature, metal- 
licity, mass, radius, and age from the IRT F/ Sp eX and 


Kec k/ HIRES spectra described in Section 2.3 


Mann 


et al. (2013b) provide spectral indices and empirical fits 


to estimate M dwarf effective temperatures from optical 
and near-IR spectra. These i ndices and fits are calib rated 
using the M dwarf sample of Boyajian et al.| (2012). Fol¬ 
lowing these methods, we estimate the .JUK-band tem¬ 
peratures of K2-21 and calculate a mean T e g- and rms 
scatter to find T e ff = 4043±375 K (±148 K systematic er¬ 
ror, ±375 K total error). This range of effective temper¬ 
atures is consistent with main-sequence d warfs sp anning 
the tr ansition between K- and M-types (Pecaut & Ma¬ 
ma jek 2013 1. We use t he met allicity calibration and solt- 
war^jof Mann et al. (2013a) to calculate K2-21’s metal¬ 
licity using the H- ana K-bands. The final metallicity is 
the mean of the H- and K-band estimates and the error 
is calculated from the quadrature sum of the measure- 



5-1.0-0.5 0.0 0.5 1.0 1.5 2.0 

Aa (arcsec) 


0.5 1.0 1.5 2.0 

Radial Distance (arcsec) 

Fig. 5.— NIRC2 K-band image and contrast curve. No stars 
with contrasts A K < 3.7 are detected with separations >0.1 arcsec 
and A K < 7.5 with separations >0.5 arcsec. 

ment error and systematic error in each band. We find 
[Fe/H = —0.11 ± 0.13 dex, approximately solar. 


Mann 


et al.|( 2013b I also provide empirical calibrations ancfsoft- 


warF’lto calculate the radii, masses, and luminosities of 
M dwarfs as a function of T e ff. Using the most con¬ 
servative T e ff errors, we calculate i?* = 0.60 ± 0.10 Rq, 
M* = 0.64±0.11 M q , and L* = 0.086±0.064 L & . These 
values, and the other fundamental parameters of the star, 
are tabulated in Tableland are used for our calculations 
of the individual planet properties. 

We also use our SpeX spectrum to estimate the spec¬ 
tral type of K2-21. We calculate the spectral type from 
several temperature sensitive molecular indices across the 
spectrum. In the optical, th e Ti05 and CaH3 indices 


(Reid et al. 1995 Gizis 19971 are covered by th e spec¬ 


trum and are calibrated lor the earliest M dwarfs (Lcpine 


et ai. 2 003 1. We follow the most recent sp ectral type cal¬ 
ibrations of these indices by Lepine et al. (2013) and esti¬ 
mate a spec tral type of K7.5±0.5. In the JN1K, the H 2 O- 


K2 index of Rojas-Ayala et al. (2012) provides a spectral 
type of M0.0±0.5. We combine the optical and infrared 
types and adopt a final spectral type of M0.0±0.5V. Our 
fundamental parameter and spectral type analyses yield 
consistent results and are also consistent with our visual 
comparison of K2-21 to M dwarf standards. We compute 
a distance modulus, p = 4.05 ± 0.19, to K2-21 by com- 
parin g the obse rve d K- band magnitude to tabulated Mk 
from 
of 65 


Pecaut & Mamajek (20131, which gives a distance 
± 6 pc. 

Our HIRES spectrum provides access to the Ha line at 
6563 A. Lines in the hydrogen Balmer series are associ¬ 
ated with magnetic activity in late-type stars and emis^ 
siondn the Ha lin e is used as a coarse age indicator (West 
et al.||2004[ |2008J) . The HIRES spectrum exhibits Ha in 
absorption, consistent with an inactive star. We also 
investigate magnetic activity in K2 -21 b y an alyz ing its 
UV emission measured by GALEX (Martin et al.||2005 ). 
The star is a weak near-UV (NUV) emitter and was only 
marginally detected (2.5cx) in the far-UV (FUV). Its low 


https://github.com/awmann/metal 


https://github.com/awmann/Teff_rad_mass_lum 
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NUV flux and marginal detection in the FUV are con¬ 
sistent with quiescent emission, similar to other nearby 


field M dwarfs (Shkolnik et al. 


Ha absorption and small U V flux are indicative of weak 


20111. For an MO dwarf, 


chromospheric activity and imply an age >1 Gyr (West 
etaL|2008 l. 


3.2. Light Curve Fitting 

We analyze the tra nsit light curves using th e same ap¬ 
proach described by Crossfield et al.| (2015|). In brief, 
we fit each planet’s transit separately using a minimiza¬ 
tion and Ma rko v Chain Mon te Carlo (MCMC) analysis 


(Foreman-Mackey et al.|2013), using JKTEBOP (|South- 
worth||201l| to model the light curves. 


When modeling the transit photometry, we adopt a 
linear limb-darkening law. While more complex (e.g. 
quadratic) limb-darkening prescriptions yield better fits 
to high SNR light curves, they are overkill given the shal¬ 
low depths of the K2-21 transits. We impose Gaussian 
priors in our analysis. For the limb-darkening parameter, 
it, we assume a distribution with center 0.56 and disper- 
sion 0.13 by referring to the values tabulated by Claret 
et al. (2012). All of the MCMC parameters show uni- 
modal distributions, and the inferred p arameters are con¬ 
sist ent with planets on circular orbits (Dawson & John- 
20121 . 


son 


Figure [2] shows the resulting photometry and best-fit 
models, and Table [2] summarizes the final values and 
uncertainties. Of particular interest are the small sizes 
and low equilibrium temperatures of K2-21b and c. K2- 
21b and c have Rp = 1.59 ± 0.43 i?® and Rp = 1.92 
± 0.53 i?®, respectively. We estimate equilibrium tem¬ 
perature assuming zero albedo according to the following 
formula: 

1 

Q. 

*~'inc 

4 (7 


T = 

eq 


where a is the Stefan-Boltzmann constant and Si nc is 
the incident stellar flux received by the planet. Plan¬ 
ets b and c have equilibrium temperatures of 5101 /q° K 
and 430±gg 0 K, respectively. These planets occupy a do¬ 
main of planet size and incident stellar radiation which 
is largely devoid of RV mass measurements and trans¬ 
mission spectra. In Section [4] we place the K2-21 system 
in the context of the current sample of small transiting 
planets. 

The transit profile constrains mean stellar density if 
one assumes a circular orbit. Since we fit each planet 
separately, we obtain two independent measurements for 
p*,cir C5 2.7^6 g cm -3 and 5. 4+35 g cm -3 , which are con¬ 
sistent with p* = 4.15 ± 0.61 g cm -3 , derived from spec¬ 
troscopy. The transit-derived stellar densities are also 
consistent with each other at the l-cr level, as expected 
for planets transiting the same star. 


3.3. Ruling Out False Positives 

Thus far, we have explained the transits seen in K2-21 
photometry with two planets around a single M0 dwarf. 
Here, we assess the possibility that other “false positive” 
interpretations can explain the photometry. One such 
explanation is that the transits are due to background 
eclipsing binaries where the eclipses are diluted by K2- 
21. Since an M6 dwarf is ss 1000 x fainter than an M0 


dwarf in the Kepler band (A My = 16.6 — 9.2 = 7.4) j^] a 
pair of maximally eclipsing M6 dwarfs would be diluted 
to 500 ppm combined with the light of K2-21. Thus, we 
consider false positive scenarios involving M0-M6 com¬ 
panions. 

Experience from the Kepler prime mission has shown 
that most multi-transiting systems (multis) are bona fide 
planets. Binary star systems have uniform inclination 
distributions and are uniformly distributed on the sky. 
Thus, the probability that two are inclined as seen from 
Earth and fall in the same software aperture is small. 


Lissauer et al. (20121 and Rowe et al. (2014) quantified 
the false positive probability tor Kepler prime multis as 
< 1%. This probability represents an upper limit for 
Campaign 3 multis, given that the field was out of the 
plane of the Galaxy (b = —57°), whereas the original 
Kepler field was near the plane of the Galaxy, where the 
density of background stars is higher. 

Despite the low probability of the transits being due 
to background sources, we searched for companion stars 
falling within the software aperture shown in Figure [4] 
No such companions separated by more than 5 arcsec 
are visible in the archival DSS images of K2-21. Our 
NIRC2 image (Figure [5) rules out companions brighter 
than AK = 4 outside of 0.15 arcsec. Over the 61 years 
between 1954 and 2015, K2-21 moved by ^5 arcsec rela¬ 
tive to distant background stars. Any background stars 
having AV > 3 within 0.15 arcsec of the current position 
of K2-21 would be detectable in the POSS1 plates. The 
low statistical probability of eclipsing binaries mimick¬ 
ing a multi-transiting system, combined with the non¬ 
detection of any background object within our software 
aperture strongly suggests that the observed transits are 
not due to background eclipsing binaries. 

We now consider the possibility that K2-21 has a 
bound companion with its own set of transiting com¬ 
panions. An M6 dwarf is 4 magnitudes fainter than an 
M0 in K-band. Consulting the contrast curve shown in 
Figure [5] we rule out companions M6 and earlier sepa¬ 
rated by more than 0.15 arcsec or 10 AU at a distance 
of 65 ± 6 pc. 

We searched for close companions using our HIRES 
spect rum of K2-21. Adopting the methodology in |Kolbl| 


et al. (20151, we searched for spectroscopic binaries in 
our HIRES spectrum. We detect no secondary set of 
lines from a star having AV < 4 mag shifted by more 
than 15 km/s relative to the lines of the primary star. 
Shifts of Av > 15 km/s correspond to orbital separations 
of < 2 AU. Thus, the HIRES spectrum rules out bound 
dwarfs M4 and earlier within 2 AU. 

While we have ruled out most of the parameter space 
where a companion star could be lurking, our search was 
not exhaustive. However, many of the remaining sce¬ 
narios involving companions with intermediate separa¬ 
tions do not pass stability considerations. A convenient 
length scale for considering dynamical interactions be¬ 
tween planets is the mutual Hill radius, 


Rh = 


Min + M out 
3 Ml 


1/3 


&Ollt 


16 Stellar param eters used in this section are drawn from|Pecaut| 
|&; Mamajek| ( |2013[ ) 
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TABLE 2 

Planet Parameters 


Parameter 

Units 

b 

c 

To 

BJD tdb - 2454833 

2156.4239+;°°“ 

2155.4708+'°[+ 

P 

d 

9 32414 +0 ' 00059 
00063 

15 5012f)~*~ 9 '00093 
lo.ouizu— 0 . 0009g 

i 

deg 

88.3±^;f 

89.08+'°° 

Rp/R * 

% 

2 60”^°' 12 
z.DU—o 15 

Q I c+0.22 

o, - lo -0.17 

T14 

hr 

2-305++ 8 

2 296 -1-0 ' 137 

T23 

hr 

2-08++ 

2 054+0- 099 

Z.UO4_ 0 174 

R*/a 

- 

0 Q43+ 0 015 

u.u^_ 0 011 

0 0244"*" 9 0100 

u.uz44:_ 0 Q050 

b 

- 

0 7Q+0.15 
u -' u —0.45 

U.O°_o 40 

u 

- 

0 45+ 013 

0.51+1° 

P*, circ 

g cm -3 

9 7+3.7 
z -'-1.6 

5.4+'j 

a 

AU 

0 0731 +0 0057 
U.LWO±_ 0 oo67 

0 1026”*" 9 0079 
U.1UZD_ 0 0094 

Rp 

R® 

1 50+0-42 

J-.oy_0.44 

L92++ 

'S'inc 

5® 

11 o +101 

5-6+1 

T oq 

K 

510+o° 

O 

O O 
tH CD 
+ 1 
O 
CO 


where M and a denote mass and semi-major axis, respec¬ 
tively. The subscripts “in” and “out” correspond to the 
inner and outer planets, respectively. The separation be¬ 
tween planet may be expressed in terms of their mutual 
Hill radius, 

A — (flout din) /77//. 


Gladman (1993) showed that two planets on initially cir¬ 
cular orbits are unstable if A < 2-\/3 ss 3.5. 

Consider the following false positive scenario: K2-21 
has a bound M5 companion at a = 5 AU that hosts 
its own two-planet system. An M5 dwarf is 5.7 magni¬ 
tudes fainter than a MO in V-band. In order to repro¬ 
duce the observed transit depths, the transits across the 
M5 dwarf must be 190x deeper to account for dilution. 
The implied planet radii would much larger: 11 R ® and 
14 i ?0 for planets b and c, respectively. In addition to 


imply planet masses of > 100 M® and A < 3.37. Such 
a system would not be stable over Gyr timescales. We 
conclude that the interpretation that K2-21 is a single 
star orbited by two planets is much more likely than in¬ 
terpretations involving companions that have evaded our 
follow up imaging and spectroscopy. 

4. DISCUSSION 

4.1. Planet Masses and Radii at Low Stellar Irradiation 

We report two small planets with sizes 1.59 ± 0.43 I?® 
and 1.92 ± 0.53 i?® orbiting K2-21, a nearby M dwarf. 
Figure [ 6 ] shows a density-radius diagram of all planets 
with 2 -<t mass and radius measurements^ Planets tran¬ 
sition from rocky compositions to possessing a signifi- 


planets are larger than 1.5 i?®, their 1-er confidence in¬ 
tervals straddle 1.5 i?® and offer a valuable probe of the 
densities and bulk compositions near the transition. 

We can estimate planet masses using mass-radius stud¬ 
ies from the prime Kepler mission. Adopting the em¬ 
pirical averag e mass-radius relationship from |Weiss fe] 


than 3 i?® 1 

Dressing & Charbonneau 

2013b 1 , this hypo- 

thetical system is dynamically unstab 

le. Such large radii 


Marcy (20141, we estimate K2-21b and c have masses of 
Mp « 4 1 W® and Mp w 5 M®, respectively. We empha¬ 
size that these mass estimates are crude given the large 
uncertanties in planet size and the large observed scat¬ 
ter in exoplanet densities for planets having Rp = 1.5— 
2.0 i?®. Assuming circular orbits, these planets would 
have RV semi-amplitudes K = 1.7 m s -1 and 2.0 m 
s -1 , respectively. Given the brightness of the host star 
(Kp = 12.3), measuring planet masses with ground-based 
spectrometers like HIRES and HARPS is possible, albeit 
challenging, assuming that th e sta r h as low levels o f jit- 
ter. As a point o f reference, Howard et al. (2013|) and 
Pepe et al.| (20131 measured the mass of Kepler-78b (Kp 


11.6), where K = 1.66 ± 0.41 m s” 

In addition to their small size, the K2-21 planets are 
also noteworthy for their relatively low levels of incident 
stellar flux. The radii of close-in planets are likely to 
be actively sculpted by atmospheric mass loss. More 
than 80% of known planets with Rp < 2 i?® receive 
> 100 times more incident stellar radiation than Earth 
(i.e. S'inc > 100 5®). Most of the exceptions orbit faint 
stars for which RV mass measurements and atmospheric 
characterization is not feasible. K2-21 is the brightest 
star known to host a planet R P < 2.0 i?®, ,Sj nc < 20 S'®, 
and P < 10 days. It joins Kepler-138 and K2-3 on the 
list of bright (J < 10.5 mag) stars known to host multi¬ 
ple small planets with low incident flux (Si nc < 20 S®). 
The K2-21 system is a laboratory to study the radii of 
planets where mass loss plays a weaker role. 


ceed 1.5 ( 

Marcy et al. 

2014 

Weiss & Marcy 2014: 

TTV measuremen 

ts of 

Rogers 

2014). 

While the most likely radii ot the K2-21 

JuiiLof-IIullei eL al. 
these nlanets (cl) is 

( 2 Ula) 
of su 


17 Data are drawn from the NASA Exoplanet Archive, July 28, 
2015, http://exoplanetarchive.ipac.caltech.edu 


ing the trend of otherwise rocky planets smaller than 
1.5 77®. It is unclear whether the low measured density 
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Radius [Earth Radii] 

Fig. 6. — Density-radius diagram for all planets with 2-cr mass 
and radius measurements, adopted from the NASA Exoplanet 
Archive. Planets with incident flux S'inc < 20 S'© (magenta tri¬ 
angles), including Kepler-138b, c, and d (red diamonds), are sys¬ 
tematically less dense than those with S'inc > 20 S© (gray circles). 
Solar system planets are indicated as cyan squares. The green 
curves indicate expected planet mass-radius curves f or pure iron , 
water, and rock compositions according to models by Zeng & Sas- 
|selov|(|2013|. Below ~1.5 R©, planets are consistent with primarily 
rocky compositions. A transition occurs at 1.5—2.0 R© and den¬ 
sity begins to decrease with radius, which has been int e rpreted as 
the onset of substantial atmospheri c accretion (|Marcy et al.|2014| 

TTV and JKV mass measure- 


|Weiss & Marcy|2014| 

ments of K2-2lb anc 


will help populate the 1.5-2.0 R© tran¬ 
sition region with low S'inc planets. Such measurements will help 
test whether low-density outliers like Kepler-138d are the result of 
limited photoevaporation, distinct formation histories of compact 
multiplanet systems, or systematic offsets in masses measured with 
TTVs vs. RVs. 

of Kepler-138 d is a product of limited photoevapora¬ 
tion, the result of different formation histories of compact 
multi-planet systems, or a systematic und erestimatio n of 


2014). 
of 


plane t masses derived by TTVs (see, e.g., Weiss & Marcy 
Combined TTV, RV, and atmospheric follow-up 
tKe K2-21 system has the potential to determine if 
small, cool, compact planets are typically low-density, 
informing models of planet formation and evolution. 

4.2. Transit Timing Variations 

Transit timing variations (TTVs) provide another av¬ 
enue by which to constrain planet properties. Measur¬ 
ing TTVs with K2 is more challenging compared to the 
prime Kepler mission because the time baseline over 
which TTVs can accumulate is much shorter («65 days 
as op posed to «4 years). Despite these challenges, Arm¬ 


strong et al. (2015) measured TTVs for K2-19b and c, 


using photometry from K2 and the ground-based NITES 
telescope. 

For the K2-21 system, the ratio of the mean orbital 
periods is P c /Pb = 1.6624, close to a 5:3 period com- 
mensurability. This proximity indicates that these plan¬ 
ets, though low in mass, could interact strongly enough 
gravitationally to produce observable TTVs. Although a 
formula for the TTVs of planets near a second order res¬ 
onance has not yet been derived in the literature, we can 
use our understanding of TTVs for pairs of planets near 
first order mean motion resonances to hypothesize the 
important parameters. The “super-period” of the TTV 
signal, assuming the pair is not dynamically in the 5:3 
resonance, will be given by 1/(5 /P c — 3/Pb) ~ 1,230 


days. The 65-day K2 time series represents only 7% of 
this super-period. 

The amplitude of the TTV signal will depend linearly 
on the mass of the perturbing planet, re lative to the 
mass of the host star (Agol et al. 2005) but also on 
the eccentricities (e) and Tongifudes of pericenter (zu) 
of each planet. For first order resonances, the TTV 
amplitude and phase depend primarily on the quanti¬ 
ties eb cos Wb — e c cos zz c and sin zzb — e c sin w c and 
not on e ach of the eccentriciti es and pericenters indi¬ 
vidually (Lithwick et al. 20121. A similar result has 
been found empirically tor second order resonances by 
fitting the TTVs of Kepler-138. Kepler-138c and d are 
near the 5:3 resonance and have poorly constrained ab¬ 
solute eccentricities and pericenters, though the differ¬ 
ences eb cos rob — e r cos zo r and eb sin Wb — e c sin zz, : are 
measured well (Jontof-Hutter et al. 2 0l5| . Therefore, for 
second order resonances, these may also be the quantities 
that control the TTV amplitude. 

To estimate the amplitude of the TTVs, we assume 
nominal masses of Mp = 4 M® and 5 M®. We choose low 
but reasonable eccentricities of eb = e c = 0.015. Based 
on the above reasoning, the maximum TTV signal should 
be achieved when the pericenters are anti-aligned, while 
it should be near minimized when the pericenters are 
aligned. 

In Figure [7] we show the predicted TTVs in these 
two limiting cases. These TTVs were calculated with 
the publicly available n-planet num erical integratio n and 
transit timing code TTVFast (Deck et al. 20141. In 
the bottom panel, the TTVs are dominated by tne syn¬ 
odic “chop ping” signal associated jwith p lanet ary con¬ 
junct ions (Nesvorny & Vokrouhlicky ,2014 Deck & Agol] 
20151. This chopping pattern oscillates o n the resonant 


timescale because of an aliasing effect (Deck & Agol 
20151. 


Assuming nominal masses and eccentricities, TTVs 
range from «1 min to «1 hr, depending on the angle 
between planet pericenters. TTVs of «1 hr should be 
detectable from the ground. A detection of TTVs would 
further constrain the orbital parameters of the K2-21 
planets and compliment RV measurements. 

4.3. Formation Scenarios 

The proximity of K2-21b and c to the 5:3 resonance 
may contain clues regarding their formation. The frac¬ 
tional deviation of the current period ratio from exact 
commensurability is 3 (P c /Pb — 5/3)/5 = —0.0025. K2- 
21 b and c may have been in resonance at some point 
in the past and may still be in resonance today. At 
this deviation from exact commensurability, the K2-21 
pair would need eccentricities of approximately 0.05 to 
be locked in the resonance, according to the estimate of 
the wi dth of the resonance derived by |Veras fe Arm itage 
(20041 and given in Equation 4.46 of that work. 

An open question is whether the close-in, small plan¬ 
ets shown to be c ommon around G, K, fc M dwarfs 
form in situ (Hansen 2014; Chiang & Laughlin 20 13t o r 
migrate inward (Izidoro et al. 2014 " Schlichting 20141. 
Multi-planet systems migrating inward are thought to 
get trapped in resonance. However, the distribution of 
period ratios is largely uniform with an excess of planets 
just outside of first order resonance and a paucity just 
inside. Theories invoking inward migration have shown 
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Time (BJD-2456988) 


Fig. 7. — TTV predictions for the two planets based on plau¬ 
sible masses and orbital elements. The case where the pericenters 
are anti-aligned (top) corresponds to a maximum TTV amplitude, 
for given masses, while the aligned case (bottom) corresponds to a 
minimum TTV amplitude. The integration lasts « 2 “super peri¬ 
ods.” Eccentricities of 0.015 are assumed. The TTVs (filled circles) 
are connected in the bottom panel for clarity. 


that it is possible for pla nets t o break out of resonance 
(Goldreich & Schlichting 2014). 

Does the close proximity of the pair to the 5:3 res¬ 
onance indicate that this system must have assembled 
via convergent migration and subsequent capture into 
resonance? There are a handful of known pairs of low- 
mass planets with period ratios that are even closer to 
the nearest low-order commensurabi lity (within a devia - 
tion of 0.001, in a fract ional sense; Fabrycky et al. 2014). 


However, as noted by Fabrycky et ai.| (]2014|), even tor- 
mation processes which produce a smooth background 
period ratio distribution will contain some pairs that ran¬ 
domly fall very close to resonance. Though the proxim¬ 
ity to resonance may dictate the dynamical evolution of 
these systems, it is not necessarily a smoking gun for 
assembly via convergent migration. 

5. CONCLUSIONS 

We report two small planet near mean motion res¬ 
onance orbiting a bright M0 star observed during the 


K2 mission. This system is a convenient laboratory for 
studying the bulk compositions and atmospheric proper¬ 
ties of small planets with low-equilibrium temperatures. 
Given that Kepler is expected to observe at least 14 fields 
in total, we expect many more exciting systems to emerge 
from the K2 mission. 
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